IPPP/12/72, DCPT/12/144 



Enhanced Gamma Ray Signals in Cosmic Proton- Wimp Collisions Due to 

Hadronization 

Spencer Chang 1 , Yu Gao 1 , Michael Spannowsky 2 

1 Institute of Theoretical Science, University of Oregon, 
Eugene, OR 97403-5203, USA 
2 Institute for Particle Physics Phenomenology, 
Department of Physics, Durham University, 
DEI 3LE, United Kingdom 

In this paper, we investigate the gamma ray signal produced from dark matter collisions with 
high energy cosmic protons. Notably, we extend past results by including important hadronization 
effects. Showering and hadronization produces a high multiplicity of photons from the decays of 
hadrons, whose rate is not suppressed by the fine structure constant. Notably, proton remnants that 
do not participate in hard scattering, can produce a large rate of photons in the forward direction. 
These effects significantly enhance the photon rate and alter the energy and angular distributions 
compared to previous results which used only parton level calculations. Due to this modification, 
the gamma ray signal from the nearby active galactic nuclei Centaurus A is potentially testable in 
future Fermi-LAT and HESS measurements, for a dark matter mass and coupling consistent with 
current XENON100 bounds. 
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I. INTRODUCTION 



The fundamental nature of dark matter remains mysterious to this day. The impressive consistency with the dark 
matter paradigm on a wide range of scales is tempered by the fact that only its universal gravitational interactions 
have been observed. To go further in our understanding, obtaining solid evidence of dark matter interactions with the 
Standard Model would be of enormous value. One promising approach to test such interactions is observing cosmic 
gamma rays. Dark matter can produce gamma rays through annihilation or decay and is currently being searched for 
by experiments such as Fermi-LAT [l[ and HESS Q . 

Another dark matter gamma ray signal that is less well known is gamma rays produced in collisions of dark matter 
with high energy cosmic particles. This indirect signal was proposed initially by Ellliot and Wells Q and has been 
reinvestigated more recently. Specifically, gamma rays from collisions between cosmic ray (CR) and dark matter 
(DM) particles have been studied in regions with concentrated dark matter and high energy cosmic rayflux, near the 
center of active galactic nuclei (AGN) [H, [j| and also the central region of the Milky Way Refs. 0t1| investigated 
the parton-level radiation and showed that the scattering between dark matter and the cosmic electrons in AGNs 
can be a promising gamma ray signal. On the other hand, Q found that proton-wimp interactions lead to a less 
promising signal due to the quark's fractional charge and momentum distribution in the proton. However, this analysis 
was a parton-only calculation and neglects important effects, such as showering and hadronization. In this study, we 
demonstrate that the hadronic shower contributes a large number of photons from hadron decays. Notably, this photon 
production is enhanced relative to the parton-only calculation due to large multiplicity and the lack of suppression 
by otQED- Furthermore, these photons have a substantially altered energy spectra, with energies extending to higher 
values compared to those from the hard scattering process. The hard scattering occurs at an energy scale related 
to the dark matter mass, which can be much lower than the total incoming proton energy. The showering from 
the rest of the proton, which does not participate in the hard scattering and is likely to carry the majority of the 
incoming energy, can emit very energetic photons in the forward direction. For AGNs, this enhancement allows the 
small fraction of protons directed towards the Earth to give a significant gamma ray contribution. To summarize, we 
find that the photons of the shower are an important modification to the gamma ray signal from proton-dark matter 
scattering which greatly enhances the rate and modifies the shape of the energy spectrum. 

In our analysis we adopt a toy model of a Majorana fcrmion dark matter, that couples to the right-handed up quark 
through a heavy charged scalar. This simple model serves as a template for models where the dark matter - cosmic 
ray collisions give an important photon signal, while other signals, like dark matter annihilation, become suppressed. 
With a heavy partner to a standard model quark, the collision process can undergo an s-channel resonance, leading 
to a large scattering cross-section. This can be consistent with bounds from direct detection experiments where the 
^keV recoil energy is much lower than the resonance energy. Interestingly, as most of the photon radiation is emitted 
on resonance, the photon spectral structure is determined by the mass difference between the DM and the heavy p/e~ 
partner, instead of the mass of DM itself. This creates freedom in the signal's energy scale that differs from that from 
DM annihilation and decay cases. 

The rest of this paper is organized as follows. In Section [TT] we discuss the resonant scattering process in the toy 
model. We take into account the allowed mass range from the latest XENON100 constraints [H, to determine allowed 
signal benchmarks. Section IIIII outlines the calculation of the gamma ray spectrum. In Section IIVI we study the 
signal from collisions off cosmic rays for the AGN Centaurus A. In Section [V] we comment on the enhancement on 
the gamma rays from dark matter collisions with diffuse cosmic rays. A summary is presented in Section IVT1 Finally, 
in the Appendix [AJ we list some important formulas for diffuse cosmic ray scattering. 



II. RESONANT COSMIC RAY-DARK MATTER SCATTERING 



The rates of cosmic ray-dark matter scattering production of gamma rays are particularly interesting when there 
is an s-channel resonance enhancement (4|. For example, a heavy scalar partner to the up quark, <fr, can mediate an 
s-channel resonance as shown in the left column of Fig[T] This should be compared with the three processes in the 
right column which illustrate the leading parton-level photon emission. 

In the galaxy, the dark matter x is non-relativistic. Thus, the condition on the up-quark's energy to hit the 
resonance is (neglecting the up quark mass) 

M 2 - M 2 

M 2 = S =(p x + Pu ) 2 =M 2 + 2E u M x or E u = E res = x (1) 

where s is the square of the center-of-mass (COM) frame energy. To realize this signal, we will consider a theory with 
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FIG. 1. Resonant s-channel diagrams. The charged u-partner, <fi, can also radiate photons. The left panel is the leading 2 to 2 
collision and the right panel shows the leading \ u ~^ X u 1 processes that dominates the hard event's photon emissions. 

a Majorana fcrmion dark matter with the interaction, 

Cmt = y xPr u 0* + h - c - ( 2 ) 

where the dark matter \ only couples right-handedly to the up quark via the scalar a colored partner to the up 
quark that carries the same electric charge. We choose the dark matter to be Majorana to avoid inducing large 
spin-independent scattering which would be strongly excluded by direct detection experiments. Furthermore, the 
right-handed coupling ensures that non-relativistic annihilation \X ~> u ^ is chirality suppressed. For proton energies 
above E res , integrating out the parton distribution function (PDF) always ensures hitting the s-channel resonance. 
Due to the enhancement at the resonance, the total high energy scattering cross-section increases as ~ y 2 . As we will 
show, the direct detection bounds are more stringent at larger coupling, which can be avoided by taking a larger mass 
gap between x and 4>- However, since this mass gap also determines the gamma ray spectra, there is a complicated 
interplay between satisfying direct detection limits and enhancing the gamma ray signal. Thus, to be concrete, we 
set y = 1 throughout this paper and will choose the mass gap to be consistent with direct detection limits. 
The 4> decay width is 

y 2 ( m 2 \ 2 

which is less than O(10~ 2 M ( f > ) in our analysis. For such a narrow width, the s-channel resonance dominates when 
kincmatically allowed and the total a xp shoots up for E p > E res . While a xp continues to grow with proton energy, 
the cosmic proton flux normally decreases as a power-law spectrum. Thus the total gamma ray signal depends on the 
energy where a xp turns up, which is determined by the mass gap between x, 4>- 

This mass gap has a lower bound from direct detection experiments since the scattering rate is enhanced in the 
squeezed limit |7|. In Fig. [5J we show the minimal mass differences allowed by recent XENON100 results [8[ at 90% 
confidence level. We have included the spin-dependent (SD) scattering and also the spin-independent (SI) scattering 
induced by the twist-2 operator Q , which are comparable in importance near the bound. Due to the resultant change 
in the energy spectrum, we cannot use XENONlOO's limit which is based on a profile likelihood. To construct our 
own limit, we use XENONlOO's hard discrimination cut acceptance shown in their Fig. 1 0] and require less than 5.3 
expected signal events, which is the 90% CL limit given their two observed events. As a cross check, our derived limit 
on the SI cross section <j x n is slightly weaker than their observed profile likelihood limit, but consistent within the 
ler expected sensitivity band in their Fig. 3 @. We refer to @, E3l f° r further details of the SD and SI cross section 
calculations. 
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As two benchmarks, we use points A and B at (m x , m$) =(300, 405) GeV and (1, 1.04) TeV, While a narrower mass 
gap can be allowed for a lower y, the scattering cross-section decreases faster than the gain from a lower resonance 
CR energy (for the assumed E~ 2 spectrum in this paper). Thus the gamma ray signal turns out to be less favorable 
with smaller couplings. At the sample points, non-relativistic annihilation XX u ^ is chirality suppressed by the 
small u-quark mass. The leading annihilation process is XX ~^ u ug, with a sub-picobarn (va) that is allowed by 
PAMELA [ll| results on the local p/p ratio. 

Although a light dark matter at multi-GeV mass may evade the direct detection bounds [la ], in our SUSY-inspired 
toy model the charged scalar partner (f> is required to have a very small mass difference from the dark matter candidate, 
in order avoid detection at LEP and LHC. While such a scenario is possible, it is beyond the scope of this paper to 
fully study the collider bound at a narrow corner of our toy model. 




200 400 600 800 1000 
M x [GeV] 



FIG. 2. Dark matter and u-partner masses allowed by 2012 XENON100 results at 90% C.L. Sample mass points are marked 
above the y — 1 bound. The XENON100 constraint with y = 0.5 (dotted) is also shown as the dotted curve for comparison. 



III. PROMPT GAMMA SPECTRUM 

When cosmic rays collide off dark matter, photons are emitted either directly from the hard scattering process 
or produced during the shower. The parton level spectrum has been studied in detail in Ref. Q. The high energy 
photons are emitted in different ways during \i e ~ ancl X>P collisions, as described below. 

In the case of x, e~ collision, the leading photon emission is through the \e~ — > x e ~l process. When the electrons 
energy is above E res , the initial state radiation (ISR) diagram takes over, with a hard ISR photon in the collinear 
direction that puts both the internal e~ and (f> propagators on resonance (under the approximation where the electron 
is masslcss). However, the total cross-section is suppressed at large incoming electron energy, where a ~ E_. As 
a result, the integrated photon spectrum over a power- law spectrum for incoming CR electrons falls sharply for 
E 1 > Mfj, - M x . 

For the case of wimp-proton collisions, the up quark's PDF in the proton takes the role of reducing the center-of- 
mass frame energy between the parton and dark matter without the necessity of extra radiation. This lifts the Eqb, 
suppression on scattering cross-section at large proton energy. Final state photons originate in the hadronic showers 
around the scattered quark, but more importantly, from the proton remnants that carry off most of the incoming 
energy. However, the PDF preference on relatively low parton momentum fraction suppresses a xp for E p ~ E res and 
below. The a xp shows a steep up-turn near E res and continues to increase with E p , which can be seen in the left panel 
of Fig. [5] Thus there is a larger contribution of high energy photons compared to the case of dark matter scattering 
off electrons. 

The prompt photons in x,P collision fall into two major categories: 

1. Final state radiation (non-remnant FSR). Photons emitted by the hard-scattered u-quark and its shower belong 
to this category. As the proton PDF ensures the s-channel <f> resonance, the final state (non-remnant) energy 
add up to the 4> mass, and the FSR photons typically has energy below 5M = — M x . However, since </> and 
X are comparable in mass, the resonance u — x system has small boost and the FSR photons can point at any 
direction and dominate the gamma ray signal at large scattering angles. The shower photons, resulting from 
hadron decays inside the hadronized jet, have in addition an enhancement due to large multiplicity and the lack 
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of suppression by the QED fine structure constant. 



2. Shower from proton remnants. Similar to the ISR photon in e ~ collisions, at Ecu > E res the remnants 
are emitted along the proton's incoming direction. These photons can be emitted at energy much higher than 
— M x but are mainly confined in this forward region. Notably, the number of the photons emitted in proton 
collisions exceed that in x, e~ ISR. At lower incoming energy Ecr rNJ E res or less, the collision does not hit 
resonance. Here, we find that the photons from remnants also are less energetic than M$ — M x but have greater 
freedom in their direction. 

We calculated the photon spectrum from x,P collisions with the Monte Carlo generator Sherpa [Hj], and its Amegic [lH 
and Photons [1 5 | p ackages for showering and photon radiation. The dark matter Lagrangian is implemented with 
the FeynRules [16( package. 



We use the hadronization model included in Ahadic [171 ] as implemented in Sherpa. It is designed to allow the study 
of deep-inelastic scattering (DIS) processefl[l8| . The resulting partons of our DIS-like hard process, px — > X + A, and 
the proton remnants are first showered and then transformed into primordial hadrons during the cluster hadronization 
process. The subsequent decays of unstable hadrons are also handled by Sherpa. In fact, most of the photons produced 
in the simulation emerge from the decay of neutral mesons, such as ir° or through final state radiation in hadron decays. 
Admittedly, the modeling for remnants is not experimentally tested for proton collisions with an neutral exotic particle. 
Thus, we warn that the theoretical uncertainties on the photon production can be large, especially in the forward 
region. 




FIG. 3. Photon spectrum at a near-threshold (125 Gev, solid) and high (10 5 GeV, dashed) proton energy. The spectra 

for parton level \ u X u l are shown in blue color, where ^ is normalized to one photon in each collision. Black curves are 
for fully showered \P with remnants and have more photons per event. For all spectra, the dark matter mass is 300 

GeV and the up quark partner is at 405 GeV. 



Fig. [3] shows the photon spectrum from one %,p collision event, comparing the fully showered x u ~^ X u case with 
the hard photon spectrum from a parton level calculation of x u ~* X u l- For the parton level x u ~^ X u l calculation 
we imposed these cuts: photon E, Pt > 1 GeV and the invariant mass between u, 7 greater than 1 GeV. Since the 
photon bremsstrahlung has a logarithmic dependence on the charged particle (u-quark) mass, these cuts may cause 
a factor of order 0(1) to the normalization of parton level radiation, and do not qualitatively impact the results. 
Here we assume an isotropic distribution of incoming protons, at 125 GeV that is close to E res (solid curves). This 
produces photons mostly below the mass gap — M x . This should be contrasted with the high energy case of 
E p = 10 5 GeV (dashed curves. Note that the gap between the fully showered and parton level spectra widens for 
higher energy incoming protons, due to extra hard photons from remnants. For E p above resonance energy, while the 
FSR photons remain below Af^ — M x , ISR and proton remnants photons continue to higher energies. 



1 We thank Stefan Hoeche for enabling Wimps as initial-state particles in the event generation in Sherpa. 
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FIG. 4. Angular distribution of photons with E-, > 1 GeV at cosmic ray energy near (left panel) and above (right panel) the 
resonance energy. 6 denotes the photon's scattering angle off the incoming proton direction. Showered spectra are normalized 
to the number of final state photons, N 1 (above 1 GeV). Parton level curves are normalized to one photon per event and their 
gap to the fully showered curve at large scattering angle is due to N-,. Note that above E re3 the ISR dominates the x e ~ ~* X e_ 7 
spectrum and emission at large angle diminishes. Furthermore, the behavior at small 9 show how proton remnants provide 
extra forward photons. For incoming particles near resonance energy (left), photons from hard event are less peaked in the 
forward direction. The x e ~ collision results are also shown for comparison. For all spectra, the dark matter mass is 300 GeV 
and the up quark (electron) partner is at 405 GeV. The parton level spectra have kinematic cuts that require both photon Pt 
and w/e - ,7 invariant mass above 1 GeV. 



In Fig. |4l we show the distribution of photons {E 1 > 1 GeV) over the scattering angle. Here the scattering angle 
8 is between the photon and the incoming proton's direction. Note that the remnants give a more pronounced peak 
in the forward direction for XiV collisions, as exemplified by the 1 TeV incoming proton. Although \e~ collisions for 
high incoming electron energy also favors a forward ISR photon, the distribution of e~ signal is dominated by <J xe - 
near E res , as seen in left panel of Fig. [5j 




£,, (0 [GeV] E y [GeV] 

FIG. 5. (left:) Dark matter - cosmic p/e~ scattering cross-sections. The parton-level \u —¥ \wy cross-section depends on 
the kinematic cut on the final state photon, besides a generic olqed suppression in comparison to that of the leading order 
\u — > x u - (right:) Unit-volume gamma ray spectra from full proton shower, parton level x u ~~ X M 7 an d X e ~ ~~ ► X e ~7- in the 
right panel, the integrated flux I(E^) is given in Eg. IA5I with a reference energy Eo = 10 GeV. The CR flux is taken as Eg. 1121 
with normalizations k p = 1 and k e ~ = 0.01. In both panels the o" xe - is scaled up by a factor of ~3 to compensate for the 
kinematic cuts. For all spectra, the dark matter mass is 300 GeV and the up quark (electron) partner is at 405 GeV. 

The cross-sections for XiP scattering are shown in Fig. [5] Scattering with an electron is also plotted as a comparison 
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for the high energy behavior of the cross section. The x( u / e ~~) ~* x( u / e ~)l processes use the same set of cuts as 
in Fig. [3J In the right panel of Fig. [5l the photon spectra from dark matter collisions with galactic cosmic rays is 
generated by convolving with a E~ 2J proton spectrum (i?~ 3 for e~) that is typical for the cosmic rays in the central 
region of our Galaxy. The normalization on the electron flux is of that of protons; see Appendix [A] for the 

definition of 1(E). Due to the lack of peaking near E res , XiP~ scattering yields a harder gamma ray signal above 
Mfj, — M x . Compared to the parton level x u ~ * X u 1 prediction, we see that the showers give a roughly two order 
of magnitude enhancement in photons for proton-dark matter scattering. As previously mentioned, this difference is 
due to the olqed suppression of the parton level cross section and a higher final state photon multiplicity N 7 from 
full hadronization/ showering. 



IV. THE AGN CASE: CENTAURUS A 



An interesting place to look for \p collision is at the center of nearby AGNs, where dark matter halo is assumed to 
exist and luminous jets provide high energy protons. We take Centaurus A for our calculations. Although the proton 
composition of the AGN jet has large modeling uncertainty, protons can make up a majority of jet particles and the 
energy output in protons, L p , from Cen. A can be more than 10 times higher than in leptons [l|. In this section, 
we compute the gamma ray signal arising from proton-dark matter collisions, with a focus on the contribution from 
proton remnants. 

We make an assumption that protons are isotropic in the AGN's 'blob' frame, similar to the jet electrons. The 
spectrum of protons are not well known. We assume that the protons also undergo Fermi acceleration and their 
isotropic spectrum in the 'blob' frame is a power-law E~ s , where the index s = 2: 



dN p 




dEdfl 



(4) 



The tilde ~ denotes variables in the boosted 'blob' frame. K p is the normalization that is determined by AGN's proton 
output. The 'blob' frame moves at a Lorentz boost factor Tb relative to the central black hole. Following the analyses 
in [j, we take Tb = 3. Boosting back to the black hole frame where the wimps are non-relativistic, the proton 
spectrum is 

dN p dN p 1 



dEdn dEdtl T B ( 1 - Pb cos 0) 
where the (3b = 0.94 is the 'blob' frame's velocity. The energy and zenith angle before/after the boost are related by 

cose= , COS0 ~ l3B a E = E-r B (l~f3 B cos6). (6) 

1 — PB COSU 

with 9 = along the jet axis. Combining Eq. [3]and|ni the black- hole frame proton spectrum is 

-1^^- ~ £ (I. V S . [F B (1 - p B cosO)]-^ , (7) 



2ndE dcos9 An \E 

which is still a power-law and has the same index as that in the 'blob' frame, while its intensity now varies with 
direction. For relativistic particles, the reference energy Eg is irrelevant and can be absorbed into the normalization. 
Keeping a nonzero mass leads to a 0(^~ 2 ) correction to the formulae above, where j p = E/m p . Integrating Eq. [7] 
with proton energy gives the total proton energy output as 



L , 



p El T B HE max I ^-^Vram) 

As Eq. [7] shows, while the flux along the jet direction is greatly enhanced, the Jacobian suppression at large angles 
oc T^ 3 . The Cen. A jet is 68° (cosf = 0.37) from the Earth; in this direction, the Jacobian from the Lorentz boost 
suppresses the proton flux by a factor of 0.14 compared to the unboostcd flux. Most protons are along the jet axis, thus 
their gamma ray contribution towards the Earth is through large angle scattering. In terms of phase-space, for these 
'along-axis' protons, a 47r integration of the Jacobian J \ J\ dfi is of order O(10 2 ) times favored by the Lorentz boost, 
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in comparison to the Jacobian integrated around an angular window Acos# < 0.1 centered on the protons pointed 
towards the Earth. Given the significant enhancement from proton remnants, the 0(10 — 10 2 ) higher photon flux can 
negate/overcome this suppression at large incoming proton energy. While the low E 1 spectrum is still dominanted by 
radiation from the protons along the jet axis, protons that point near to the Earth also make considerable contribution 
to the gamma ray signal, especially when E 1 > — M x . 

Admittedly the Lorentz boosted 'blob frame' is a simplistic picture for the protons inside the AGN jet. If the 
AGN jet is more collimatcd than our assumption, less protons would point towards the Earth and the photons in the 
forward region would play a less important role. To illustrate this uncertainty in the AGN proton distribution, our 
signal prediction will also be shown with just the 'along-axis' protons with cos8 p > 0.8. 

Including radiation from protons in all directions, the photon flux towards Earth is 



dE 7 dn 



1 Sdm 
R 2 M v 



dilp J dE p <t{E p , 



dN n 



dN~, 



dEpdilp 



d^J^y di^\^y SC. 



(10) 



Throughout this paper we denote angular- integrated cosmic ray flux as <f) and its angular differential form as dcf)/dtt. 
The direction of the incoming proton 9 p = {9 p , 4>p} is not limited to the vicinity of the jet axis. 6 1SC denotes the 'real' 
photon scattering angle in a frame where the proton momentum is along the z-axis. 7iSC is determined by the proton 
direction 9 p and the Earth's direction 6®. R = 3.7 Mpc is the Earth's distance to Cen. A. The integrand in Eq. [TOl 
determines the contribution from protons at angle 6 p off the AGN jet axis. Sdm — (Px( r ) ' r ) l& the dark matter halo 



density integrated over the distance range where collisions occur; 5dm at Cen. A can be as high as IO^Mq/pc 2 
For the jet output in protons, we use L p =lxl0 45 erg s _1 and an energy range [E m i n , E max ] = [10, 10 7 ] GeV 



dN 



denotes the final state photon distribution from an average collision event. 



dN 
dE-,d{1.-, _ sc 



4] 



is normalized to the 



dE^d^l-y 

total number of photons above 1 GeV per collision and is Monte Carlo generated. 

As shown in Fig. [5J the protons along the jet axis (cos 6 P <~ 1) suffices for parton level radiation (shown in blue dotted 
contours). The showered spectra (shown in black solid contours) shows that the forward photons, with cos9 p ~ 0.37, 
mostly from proton remnants, account for a significant portion or even the majority of the signal, especially at large 
E~(. Even at low E 11 contributions from protons along the jet axis are still significant. The resulting gamma ray 
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FIG. 6. Emission intensity contours in terms of proton energy and proton's angle off the jet axis. Black and blue contours 
denote the showered (solid) and parton level (dotted) calculations, respectively. The azimuthal angle <f> p — in both panels, 
i.e. protons are in the Earth - jet axis plane. The Earth's direction is cos(9 p =0.37. Emission intensity / is the integrand in 
Eg. I lOt in units of pb-GeV~ 2 -sr~ 2 • K p \ Eo=10 GcV . For all spectra, the dark matter mass is 300 GeV and the up quark partner 
is at 405 GeV. 

signal at the Earth is plotted in Fig. [7] in black. The signal is noticeably enhanced over the parton-level calculation 
shown in blue. Furthermore, the shape is substantially altered. At low energies, E 1 < — M x , the spectra is softer 
than the parton-level result and has less of a peaking structure. At high energies, E~ ( > M$ ~ M x , the fully showered 
spectrum is a power law at high energy which receives a significant contribution from protons along the jet axis. In 
comparison, the parton level photons drop abruptly after reaching — M x . Thus, taking into hadronization and 
showering has both significantly enhanced the signal and altered its spectral shape. Note that due to numeric stability 
issues in Monte Carlo, we do not plot E 7 above 1 TcV, yet the E~ 2 power-law spectrum is expected to extend to 
higher energy, as a fragmentation from the the total incoming proton energy 
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E y [GeV] 

FIG. 7. XjP collision induced gamma ray signal from Cen. A, for sample point A (solid) and B (dotted). The full collision 
and parton level E 2 -^ spectra are plotted in black and blue colors, respectively. The orange curves show the component of the 
fully showered spectra that originate from protons along the AGN jet axis, cos 8 P > 0.8. Signal levels assume optimistic AGN 
parameters: L p — 10 4j erg s _1 and 8dm = 2 x lO 11 M p C - 2 [|. Fermi Q 

and Hess [2| measurements are shown for comparison. 



In Fig. [7] , we have also plotted the signal component from the protons that are just along the AGN jet axis (shown 
in orange), integrating over cos^j, > 0.8 in Eq. [TO] This demonstrates the situation for a more highly collimated AGN 
jet, where only large-angle scattering gamma rays contribute. Thus, the high energy tail above Ma, — M x is sensitive 
to the theoretical uncertainties of the AGN jet angular distribution as well as the photon contribution from proton 
remnants. As these curves show, for large angle scattering, the major difference is that such photons drop much more 
abruptly when Ey approaches to Ma, — M x . The parton-level curves, while their normalizations arc suppressed by 
olqed/N^, also demonstrate this high energy behavior of photons from hard scattering. 

The overall signal level scales linearly with 8dm and the AGN's energy output in protons. In illustrating the gamma 
ray signal we assumed an optimistic scenario with regards to the values of the dark matter density, the AGN's proton 
energy output and the interaction coupling. The resulting gamma ray signal level for sample point A is comparable 
to the uncertainties in the Fermi data and future observation may constrain the coupling to lower values. More 
optimistically, with further enhancements to the dark matter signal, the high energy tail could explain the HESS data 
points without modifying the lower energy Fermi points and thus resolve the discrepancy in power law observed by 
the HESS and Fermi-LAT analyses [J. 



V. DIFFUSE PROTONS 



For an isotropic distribution of protons, e.g. the diffuse protons inside the Milky Way, there is no prefered direction 
and the contribution from proton remnants are present in the 47r-averaged prompt spectrum. However, due to 
relatively low CR flux inside the Milky Way plus a high energy threshold for resonance scattering, the gamma ray 
signal is much below galactic background levels. In this section we only describe the calculations with two template 
profiles of galactic protons. 



The photon signal is given by, 



drO^J. / dE v <j YV —?- 11 



dEj J Mx R - dE ri dE- 



where p(r) is the dark matter halo density, is the cosmic ray flux. At the center of the Milky Way galaxy, these 
fluxes can be parametrized [l!| as 



-2.7 



— - — — = k v ( ^ P 1 GeV 1 cm 2 s 1 sr 1 for protons 

dE p dn p \GcVj 1 

,^ e ' - fc e - ( GeV-WW- 1 for electrons (12) 
dE e -dn \GeY J 
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Note: the flux normalizations k p / e - in lower-case are not to be confused with that of the AGN jet. Since the diffuse 
spectrum is isotropic, the forward photons are readily present and their spectrum is illustrated in Fig. [3J 

For protons, its power law index only varies slightly during propagation and the spatial and energy parts in Eq. QT] 
can be separated, 



dE^dfl 



(9) = J(9) ■ I(E 7 ) 



where J(9) integrates over the dark matter distribution along the direction 9, while I(E 7 ] 
spectrum convoluted with the proton energy spectrum. For details see Appendix 



(13) 

is the prompt gamma 



CR model 


a p 


M x ■ J(0) Central 


M x ■ J{6) Inner 


Plain diffusion [20] 


-2.68 


3.8 


0.65 


Diffusion reacc. [21] 


-2.75 


9.9 


1.5 



TABLE I. Angular averaged M x ■ J in central 



< 1° ) and inner (\l\ < 



< 8°) galactic regions. a p is the power index 



of the proton flux. J is evaluated with reference energy Eo at 10 GeV. See Eg. IA4l for definitions. The dark matter profile is 
given in Eq. [14] M x ■ J(6) values are in the unit of 10 28 s" 1 m~ 4 sr~ 1 . 



We choose two template Galprop CR profiles, the plain diffusion model 999726 [20j and diffusion-reacceleration 
model 599278 [2l[ to calculate the integrated strength of gamma ray source, [ J(9)d9 in Fermi's angular windows (22|, 
as shown in Tab. |TJ For the dark matter halo we pick the Einasto profile [23| as an example for cuspy dark matter 
distribution, 



Px = Pe* 



K( r,a - r ©)/ r ?] 



(14) 



where a = 1.7, r = 8.3, r s = 25 kpc and the local halo density pq = 0.3 GeV/cm 3 . 

Inside the Milky Way, however, the proton-dark matter scattering is at an disadvantage to due the relatively low 
level of cosmic ray flux. Small mass splitting between <f>, x may face increasingly stringent constraint from direct 
detection experiments. As Eq. Q]and Fig. [2] illustrate, only cosmic protons of O(10 — 100) GeV or above significantly 
contribute to the gamma ray signal for our toy Lagrangian. As the result the p, x collision signal is dwarfed in 
comparison to the astrophysical background. 



VI. SUMMARY 



In this paper we investigated the gamma ray signal from the collision between dark matter and high energy cosmic 
ray protons, including the effects of hadronization and showering. This extends previous parton-level only calculations, 
with a substantial modification of the rate and energy spectrum of the photons. The rate is significantly enhanced, 
since the photons produced in hadronic decays have a high multiplicity and are not suppressed by the fine structure 
constant. In particular, we emphasize the contribution from the energetic proton remnants, which boost the high 
energy tail of the gamma ray spectra. To illustrate the shower enhancement to a parton level photon radiation, we 
implemented a simple Majorana fermion dark matter that couples right-handedly to the up-quark, to avoid large 
annihilation rates. We used the latest XENON100 limits to select viable test masses for the dark matter and scalar 
u partner for the cosmic signal from Cen A as an AGN candidate, and the case of diffuse protons in our Milky Way 
galaxy. 

We use the Monte Carlo generator Shcrpa to simulate one-sided proton remnant in a x,p collision event and 
subsequent showering. Due to limited choice of generators that allow remnant showering with exotic particle beams, 
it is of interest to further test the photon radiation from proton remnants with alternative means of calculation. This 
will help to determine whether there are large theory uncertainties on this gamma ray signal. An additional source 
of uncertainty is the amount of AGN protons which are pointed towards the Earth. This also affects the high energy 
photons and thus, improvements in AGN modeling will also help pin down this part of the spectrum. 

To summarize, we find that hadronization and showering substantially enhance the signal and in particular, the 
proton remnants significantly enhance the signal rate for the energy range E 1 > M,p — M x , making the protons pointed 
directly towards the Earth the major contributor at large photon energy. The gamma ray signals from proton-dark 
matter collision is found to be at a level which could be potentially constrained by forthcoming Fermi measurements. 
In contrast, collisions with diffuse protons inside the Milky Way suffers from the low proton flux at energies high 
enough to reach resonance scattering. However, if the dark matter candidate and the mediator are light the signal 
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from diffuse protons may become more prominent. 
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Appendix A: Diffuse cosmic ray scattering 

The prompt photon flux at the Earth is an integral over collision sources inside the observation angular cone Af2, 

d(j) 7 f dV dN 7 



dE 1 J 47rr 2 dE 7 dV 



dr ±n r°° d^ P{r) d*(E P ,E 7 ) (A1) 



47r|r| 2 J Ei M x dE p " ' dE y 

where r = (r, 6) with 6 in galactic coordinates, p x is the dark matter density and is the Air averaged galactic 
proton flux. 

Only the high-energy protons are relevant in our study and the spectrum can be parametrized as a power-law, 

where the power index s grows slightly as protons lose energy during their propagation to the outer region of the 
galaxy. In the Galprop models we adopted, the variation Ss ~ 10 -2 in the power-law index is insignificant and s can 
be approximated as a constant. 

Thus the energy integral can be separated from the spatial one and Eq. IA1I can be written into 

^ = J < e >- 7 <^ < A3 > 

where .7(9) = / °° *fe^l^(E , r ,e), (A4) 

The energy integral I{E 1 ) gives the shape of the prompt gamma ray spectrum and is independent from astrophysics. 
In the case with hadronic shower, the differential cross-section in Eg. I A5I is replaced with 

da(E p ,E 7 ) _ fJp ^dN 7 



atot(E p )—t. (A6) 



dEry dE>y 
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